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ABSTRACT 

A highly convergent "relaxation" method for the 
inversion of the full radiative transfer equation 
has been developed. The results of the iterative 
solution indicate that convergence can be achieved 
over a wide range of initial "guesses," enabling 
the temperature profile of a relatively unknown 
atmosphere to be unambiguously determined. The 
method is illustrated by examples for the outgoing 
radiance in the earth's atmosphere for the region 
of the 4.3y CO 2 band, but can be similarly applied 
in other frequency ranges. 


1. INTRODUCTION 

The total upward radiance I(VjP), arriving vertically in a small solid 
angle (w s» 0) at a pressure level P of an atmosphere in a given wave-number 
interval v, is obtained by formal integration of the fundamental radiative 
transfer equation (RTE) , subject to a black-body boundary condition, and is 
conveniently reduced to read 

P 

I(v,P) = e B[v, T(P 0 )]t(v,P 0 ) + / B[v,T(P)] — ( ^ P) dP, (1) 

p 0 

where B is the Planck function, t is the transmission function dependent on 
v,P, and the temperature distribution T(P) between P and P, Pq is the pressure 
at the lower boundary, and e is the emissivity at Pq. 

Following a proposal made by Kaplan (1959) that the vertical distribution 
of atmospheric temperature, T(P) , can be inferred from measurements of outgoing 
radiance, I(vi, "P ) , i = 1, 2, ... n, a number of sophisticated but inconclusive 
attempts have been made to solve Eq. (1) for T(P). The history of these at- 
tempts is extensively discussed and referenced in a series of articles by Wark 
and Fleming (1966). A review of these attempts further discloses the need for 
discovering and employing to advantage other physical and mathematical pro- 
perties of the RTE that relate the measured radiance to the desired atmospheric 
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parameter, and to a general method of solution. 

Past attempts to solve this problem, by linearizing the Planck function 
and reducing Eq. (1) to an n x n system of linear equations, failed because the 
resulting system is always "ill-conditioned." In fact, our extensive studies 
of the above method applied to the V 3 CO 2 band show that, without introducing 
any damping functions, the iterative solution converges only if the initial 
guess is made accurate to at least 10 -8 in T(P). However, this author's past 
experience (Chahine, 1965) with the same family of Eq. (1) strongly suggests 
that the RTE in its integral form is "well-behaved," and that its reduction to 
a linear system of equations is improper. 

The purpose of this work is to explore briefly an important property of 
the radiance emitted vertically from different layers of an atmosphere, and to 
apply it to the numerical solution of the complete RTE for the determination of 
the temperature profile. 

2. GENERAL FORMULATION 


The physical properties of outgoing radiation from various layers of an 
atmosphere are such that at any pressure level P within a wide range of pressures, 
P 2 <P<Pl, the total upward radiance can be related (e.g., from the mean value 
theorem) to the Planck function at P by 

Kv,P) », B[u,T(P)] 

1 r * (2) 
IMv/F) B[v,Tl(P)] 

where T(P) and T^P) are two different temperature profiles in the same atmos- 
phere, resulting in radiance values I and I 1 , respectively. This general fea- 
ture is directly related to the fact that over a wide range of temperature 
profiles, the dependence of x(v,P) on changes in T(P) is negligible compared to 
that of the Planck function. Thus , if we select a set of i frequencies that 
satisfy Eq. (2) at an appropriate set of i pressure levels, we can readily use 
the above approximation to solve Eq. (1). In practice, it is more convenient 
to proceed first by factorizing the Planck function in terms of a function of 
temperature only and then use the latter as a universal variable of iteration. 

To factorize the temperature in the Planck function we write 


B(v,T) = av 3 /(e bv/T - 1) 
= X(T) 3 (v) 6 (v ,T) , 


(3) 

(4) 


with 6 (v,T) — 0(1.0). Since the temperature dependence of 
^ -bv/T, 


we 


squares fit 


replace 
as 


the argument -bv/T, within a small v-T 

2 

dT dv = 0, (j 


V" 

T 1 

c 

Tbv 

c 3 

i 

l 

IT‘ 

C 1 ~ C 2 V " “ 


Eq. (3) appears as 
range, by its least- 


= 1, 2, 3.) (5) 


The limxts of integration reflect roughly the frequency range in use and the 
expected temperature variations in the atmosphere. Solution of the above 
3 x 3 system of equations gives 
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ln(T"/T') v" + v' 

C 1 — D rp t t rji * 2 ’ 


(6-a) 


and 


C 2 b 


C 3 . b 


ln(T"/T') 

!Jl| ! IJI I 


v" + v' 


(6-b) 


(6-c) 


The scaled universal variable of iteration is then 


X(T) = exp 
The factorized Planck function becomes 




(7) 


51^11 = X(T) 

av 


exp 


(-^C 1+ i3) 


1 - e 


-bv/T 


( 8 ) 


in which the variation of X(T) with temperature can be several orders of magni- 
tude larger than that of the second factor, as in the 4.3y region of our example. 

By substituting Eq. (8) into Eq. (1) and grouping, the RTE takes its final 
form as 


i(v,p) = x[t(p 0 )] A(v,P 0 ) + J X [T(P)] K(v,P) dp, 


(9) 


0 


in which the relatively weak dependence of A and K on changes in T(P) is not 
explicitly shown. Thus, the problem reduces now to solving^Eq. (9) for Xi when 
P is given together with a selected set of radiance values Ii (v,!*) measured at 
i frequencies. 

3. METHOD OF SOLUTION 

We propose to solve Eq. (9) by iteration: We make_a first guess T°(P) for 
X±, and using this as input to Eq. (9), evaluate I? (v,P) for the given atmosphere. 
The mechanism for generating the iterative solution is obtained from Eqs. (8) 
and (2) as prescribed in the following relaxation equation: 


x“ - x"- 1 A, 

1 1 J.n-1 

i 


( 10 ) 


where n is the order of iteration. From the above equation we determine X^, 
then T-^CP-l), substitute it into Eq. (9), and evaluate 1^. Using Eqs. (10) and 
(7) again, we evaluate xf, then T^Pi), and repeat the same procedure for the 
following iterations. It is important to point out here that when X^(T) is 
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reduced to temperature according to Eq. (7), the resulting values are tempera- 
tures at specific altitudes and are not average values taken over an atmos- 
pheric layer. 

The rate of convergence is determined from the residuals 



i 


or their average value 



The solution converges when R^ -*■ 0 for all v.. 
4. APPLICATIONS AND DISCUSSION 


(ID 


( 12 ) 


In the remaining part of this work we will demonstrate the way the above 
scheme can be used to determine the temperature distribution in the earth's 
atmosphere from the spectral distribution of vertically outgoing radiance in 
the region of the 4.3y C0 2 band. 

4.1 Procedure 

To set up the problem, we assume a temperature profile T(P), determine 
the transmission function x , and use an accurate quadrature formula to evaluate 
the radiance I^(v,P) for e = 1, P/Pq = 0.009, and for the 39 frequencies from 
2180 (5) 2370 curl (because our transmission subroutine had been designed to 
give x(v,P) at intervals of v equal to 5 cm~l in the above range). The aim now 
is to take the radiance at a set of selected frequencies as simulated experimen— 
tal data and to investigate how accurately we can retrieve the temperature pro- 
file T(P) . 


The results presented here were obtained according to the following 
sequence: 


1. Select a set of 10 sounding frequencies ( Vi = 2180, 2235, 2245, 2260, 
2290, 2295, 2305, 2315, 2370, 2360), corresponding to 10 pressure levels 
(Pi/P()_ = !- 0 » 0.85, 0.60, 0.40, 0.20, 0.10, 0.04, 0.025, 0.0135, 0.009). The 
selection of this set is not unique and the choice of its size is optional. 

The most efficient sounding frequencies, for numerical computation purposes, 
are those that satisfy Eq. (2) in the region where their contribution to the 
integrand of Eq. (9) is substantial. 

n 2. Make an initial guess T®(P) and solve Eq. (9); check the residuals 
Ri and \R^ av , and satisfy convergence. 

3. Compare the results with the exact answer, T(P). 

Other pertinent details of the above sequence will be presented along with 
our examination of the following results. 
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,4.2 Discussion of Numerical Results 

Figure 1 shows the results of the fifth iteration using the U.S. Standard 
Atmosphere as initial guess. The temperature profile computed at 10 pressure 
levels (not layers) has an average accuracy of 1°K. In the calculations, the 
transmission function t was left unchanged from that of the Standard Atmosphere 
temperature profile. When x was recomputed after every iteration, to account 
for changes in the temperature profile, the results improved slightly. When 
other sets of were used as sounding frequencies, the corresponding results 
compared favorably with those shown in Fig. 1, indicating the general prevalence 
of the feature given in Eq. (2). The numerical evaluation of the integral in 
Eq. (9) was performed by using a modified Simpson's rule and a second-order 
interpolation formula for the intermediate values of temperature. 

The rate of convergence at each of the selected 10 sounding frequencies, 
as well as that at each of the remaining 29v, was checked after every iteration. 
The average rates of convergence ^Ri^ av , listed in Table I, were taken over all 
39v. Examination of this list clearly reveals the two most important character- 
istic features of the iterative solution: stability and convergence. 

A similar examination of the average rates of convergence of the tempera- 
ture profile T n (P) indicates that an average accuracy of 1°K is achieved when 
^R n )av < 3%, and that an average accuracy of 2-3°K can be achieved when 5% < 
<^R n ^av < 7%. The above results offer an initial indication of the effect on 
the accuracy of the final profile of errors in the measurements of radiance 
values . 

Figure 2 shows the results of the eighth iteration using the same sounding 
frequencies as in Fig. 1, but using as initial guess any one of three iso- 
thermal temperature profiles, 280, 240, and 200°K. The transmission function 
x was recomputed during the first five iterations only and was kept unchanged 
afterward. After the fifth iteration, the results for all three (very differ- 
ent) initial guesses converged to the same profile, and became indistinguish- 
able on the scale of Fig. 2. The fact that the same set of sounding frequencies 
gave equally good results for all the above cases demonstrates the strong per- 
sistence of the basic feature given in Eq. (2). 

The rates of convergence for these three initial guesses were remarkably 
similar and the residuals decreased very sharply. The results listed in Table 
I show a spectacular example of convergence for the case of T®(P) = 280°K; 
there, the average residual decreased from an initial value of 920.10% to 
21.74% in just one iteration. 

It was further noticed that convergence occurred first at the lower 
levels of the atmosphere, then spread upward. Thus, higher order iterations 
tended to improve convergence at the upper levels at the expense of accumu- 
lating truncation errors at the lower levels. We note in conclusion that the 
numerical results shown here do not represent the best in the art of computa- 
tional refinements; we believe this method is capable of far greater accuracy. 


447 



PRESSURE RATIO P/P t 


RADIOMETRIC DETERMINATION OF ATMOSPHERIC TEMPERATURE PROFILES 



TEMPERATURE T t °K 

Fig. 1. Temperature profile, from 10 sounding frequencies, at the fifth 
iteration using the U.S. Standard Atmosphere as initial guess. 

• •• > iterative solution. — • — •, exact profile. 

__ __ , initial guess. 
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TABLE I. Average Residuals R n „ . % 

av 


Initial 

guess 




Order of : 

iteration, n 




T°(P) 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

U.S. 

Standard 

Atmosphere 

26.30 

6.97 

2.25 

1.52 

1.02 

0.73 

0.54 

0.40 

* 

* 

280°K 

920.10 

21.74 

10.11 

5.49 

2.88 

1.35 

0.81 

0.54 

0.42 

* 

240°K 

80.23 

21.86 

10.33 

5.49 

2.85 

1.33 

0.79 

0.54 

0.43 

* 

200°K 

90.89 

22.42 

10.76 

5.53 

2.84 

1.31 

0.79 

0.55 

0.46 

0.41 


*The relative error caused by the numerical quadrature precludes making 
the average residuals less than 0.39%. 


5. CONCLUSION 

The great advantage of the approach discussed here lies in the fact that 
it has exploited a dominant physical property of radiance occurring over a 
wide range of frequencies, and related this property to an efficient relaxa- 
tion method of solution, and that it can be easily adapted for use in different 
frequency ranges, such as the 15 y CO 2 band and the microwave regions. 

The resulting inversion scheme provides a reliable tool for practical 
sounding of temperature with the accuracy required for use in numerical 
weather forecasting in the earth's atmosphere. Perhaps more important, this 
method is also suitable for unambiguous determination of temperature profiles 
of the comparatively unknown atmospheres of other planets. The time required 
for the reduction of a set of measured radiance values to the temperature pro- 
file shown here is very small, and the calculations are simple and can even be 
carried out on a desk calculator. Thus, a readily feasible system for the re- 
duction of the measured radiance values can be arranged to produce temperature 
profiles in "real time" as data are accumulated. 
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